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ABSTRACT 

We detected X-ray emission from the Vela-Uke pulsar B 1 800-2 1 and resolved its synchrotron nebula with the 
Chandra X-ray Observatory. The pulsar's flux is fps, = (1 .4 ± 0.2) x 10"'"* ergs cm"^ s"' in the 1-6 keV band. 
Its spectrum can be described by a two-component power law (PL) -i- blackbody model, suggesting a mixture 
of thermal and magnetospheric emission. For a plausible hydrogen column density wh = 1.4 x 10^^ cm~^, the 
PL component has a slope Fpsr = 1.4 ±0.6 and a luminosity L™"* « 4 x 10-"(ii/4kpc)^ ergs s"'. The properties 
of the thermal component (kT ^ 0.1-0.3 keV, Lp°' ~ lO^'-lO^"^ ergs s"') are very poorly constrained because 
of the strong interstellar absorption. The compact, k, 7" x 4", inner pulsar-wind nebula (PWN), elongated 
perpendicular to the pulsar's proper motion, is immersed in a fainter asymmetric emission. The observed flux 
of the PWN, including its fainter component, is Fpwn = (5.5 ±0.6) x 10"'"^ ergs cm"^ s~' in the 1-8 keV band. 
The PWN spectrum can be described by a PL model: Fpwn — 1.6 ±0.3, Lpwn ~ 1.6 x 10^^(6//4kpc)^ ergs 
s"', for hh = 1.4 x 10^^ cm"^. The elongation of the inner PWN with respect to the direction of the pulsar's 
proper motion suggests that its X-ray emission emerges from a torus associated with the termination shock 
in the equatorial pulsar wind. Such an interpretation provides additional support for the alignment between 
the pulsar's velocity and the spin axis found for several other pulsars. The asymmetry in the fainter, more 
extended emission could be attributed to nonuniform properties of the ambient medium. A lack of any signs of 
bow-shock morphology suggests that the pulsar moves subsonically in high-pressure interiors of a supernova 
remnant. However, similar to a few other Vela-like pulsars, no supernova remnant is seen in the Chandra image, 
possibly because its soft X-ray emission is absorbed by the interstellar medium. The inferred PWN-pulsar 
properties (e.g., the PWN X-ray efficiency, Lp^^/E ~ lO"'*; the luminosity ratio, Lpwn/ips"* ~ 4; the pulsar 
wind pressure at the termination shock, p, ~ 10"' ergs cm""*) are very similar to those of other subsonically 
moving Vela-like objects detected with Chandra (Lp„,/E - lO-^^-lO"^ ^ Lpw„/L™"* ^ 5, p,, - lO-'^-lO"** 
ergs cm"'). 

Subject headings: pulsars: individual (PSR B 1 800-2 1 =11803-2137) — stars: neutron — X-rays: stars 



1. INTRODUCTION 

Observations with the Chandra X-ray Observatory have 
shown that many young pulsars (r < 30 kyrs) power X-ray 
nebulae (e.g., Kaspi et al. 2004; Gaensler & Slane 2006, and 
references therein). The observed X-ray emission is produced 
by relativistic particles gyrating in the magnetic field down- 
stream of the termination shock in the pulsar wind (Ken- 
nel & Coroniti 1994; Arons 2004). The innermost parts of 
many X-ray pulsar-wind nebulae (PWNe) show axisymmet- 
ric morphologies, including toroidal structures and jets along 
the pulsar's spin axis. However, even pulsars with very simi- 
lar spin-down properties (such as the period, P, period deriva- 
tive, P, spin-down power, E = 4ttIPP~^, and spin-down age, 
r = P/2P) can produce PWNe of quite different shapes and 
sizes. These differences can occur for a number of reasons. 
For instance, upstream of the termination shock (i.e., closer to 
the pulsar) the properties of the wind should be sensitive to the 
angle between the pulsars's magnetic and spin axes and the 
wind magnetization. These factors could affect the strength 
of the termination shock and, consequently, the temperature 
and radiation efficiency of the post-shock flow. On the other 
hand, the properties of the post-shock flow should also de- 
pend on density and temperature of the ambient medium. For 
instance, these parameters may vary significantly for PWNe 
residing inside supernova remnants (SNRs) of different ages. 
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Finally, the direction of the kick acquired by the neutron star 
(NS) at birth with respect to its spin axis, and the NS velocity 
relative to the ambient medium also affect the PWN appear- 
ance. For instance, if the pulsar moves with a high velocity, 
the ram pressure caused by its motion can exceed the am- 
bient gas pressure, resulting in a bow-shock PWN, such as 
"the Mouse" PWN around PSR J1747-2958 (Gaensler et al. 
2004) and "the Duck" PWN around PSR B 1757-24 (Kaspi et 
al. 2001a). Studying X-ray bright, nearby PWNe helps to dis- 
entangle various effects, understand their impact on the PWN 
structure, and probe the properties of the ambient medium. 

The well-known nearby (d « 300 pc) Vela pulsar and its 
X-ray PWN (Pavlov et al. 200 lab, 2003; Helfand et al. 2001) 
have become an archetype for young (r 10-30 kyrs) and 
energetic (E ^ lO-'^-lO"'^ ergs s"') pulsars powering X-ray 
PWNe. These pulsars are also interesting because at this age 
emission from the hot NS surface becomes observable in X- 
rays as a thermal "hump" on top of the flat non-thermal spec- 
trum (e.g., Kargaltsev & Pavlov 2006). 

Based on its age, t = 16 kyr, and spin-down power, E = 
2.2 X lO^*" ergs s"', the radio pulsar B1800-21 (hereafter 
B1800) is similar to the Vela pulsar, but it is more distant 
(the dispersion measure distance is = 3.8 ± 0.4 kpc, accord- 
ing to Cordes & Lazio 2002). Its spin-down flux, E/Aird^ = 
1.2 X 10~'^d^^ ergs s"' cm"^, where d4 = ii/(4kpc), places it 
among the top 20 pulsars ranked by this parameter The pul- 
sar is projected near the western boundary of the radio SNR 
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G8. 7-0.1 (W30), surrounded by HII regions and molecular 
gas (Odegard 1986, and references therein). However, the as- 
sociation between B1800 and G8. 7-0.1 has been considered 
doubtful because it would require a very large velocity of the 
pulsar if it was bom near the apparent SNR center (Frail et al. 
1994). 

Being young and energetic, B1800 should be surrounded 
by a PWN, and both the pulsar and the PWN should be de- 
tectable in X-rays. Based on a 10 ks ROSAT PSPC obser- 
vation, Finley & Ogelman (1994) have reported a faint X- 
ray source (PSPC count rate of 1 .5 ± 0.5 counts ks~*) located 
near the radio pulsar position and attributed this emission to 
B1800. The low count rate and insufficient angular resolu- 
tion of ROSATPSPC did not allow Finley & Ogebnan (1994) 
to resolve a compact PWN. They, however, detected diffuse 
X-ray emission ~ 30' northeast of the pulsar and attributed it 
to the G8. 7-0.1 SNR. These authors defend the association 
between the pulsar and G8. 7-0.1 by suggesting that the SN 
explosion, which produced B 1800 and G8 .7-0.1, occurred in 
or near a molecular cloud, very close to the current position of 
the pulsar, and the SNR have been "blown out" eastward into 
a low-density interstellar medium (ISM). However, a recent 
proper motion measurement (Brisken et al. 2006) has shown 
that the pulsar was bom outside the currently seen SNR, and it 
moves more nearly toward the center of G8.7-0. 1 rather than 
away from it, which makes their association very unlikely. 
This measurement also essentially rules out the association 
between the pulsar and the newly discovered SNR candidate 
G8.3 1-0.09 (Brogan et al. 2006). 

B1800 is of particular interest because it is located in the 
vicinity of the TeV 7-ray source HESS J1804-216 (Aharo- 
nian et al. 2006), and it may supply ultrarelativistic electrons 
that generate the TeV radiation by upscattering the photons of 
the cosmic microwave background. We will discuss the possi- 
ble connection between B1800 and HESS J1804-216, as weU 
as other candidate X-ray counterparts of the TeV source, in a 
separate paper. 

In this paper, we describe the results of a Chandra observa- 
tion of PSR B 1800-21 and its synchrotron nebula. The details 
of the observation and the data analysis are presented in §2. 
We compare the X-ray properties of the B1800 pulsar and its 
PWN with those of other Vela-like pulsars-PWNe and discuss 
implications of our findings in §3. Our main results are sum- 
marized in §4. 

2. OBSERVATIONS AND DATA ANALYSIS 

We observed B1800 with the Advanced CCD Imaging 
Spectrometer (ACIS) on board Chandra on 2005 May 4. The 
useful scientific exposure time was 30,236 s. The observation 
was carried out in Faint mode, and the pulsar was imaged on 
S3 chip, « 7"5 from the aim point. The detector was oper- 
ated in Full Frame mode, which provides time resolution of 
3.2 seconds. The data were reduced using the Chandra Inter- 
active Analysis of Observations (CIAO) software (ver. 3.2.1; 
CALDB ver. 3.0.3). 

2.1. Image 

Figure 1 shows the ACIS-S3 image of the region around 
B 1 800. An extended X-ray source is clearly seen in the image 
around R.A.= 18''03'"51 ?432, decl.= -21°37'07;'45 (these are 
the coordinates of the center of the brightest pixel). The dif- 
ference of 0."3 between this postion and the radio position' 

' Note that the radio timing position of B1800, R.A.= 18''03"51?333(8), 
decl.= -21°36'27"(3) (Manchester et al. 2005), differs by » 40" from the 



of B 1800 (Brisken et al. 2006) is smaller than the error in ab- 
solute Chandra astrometry (0."6 at the 90% confidence level). 
The close match between the X-ray and radio positions and 
the extended morphology of the observed X-ray emission al- 
low us to conclude that we detected the X-ray emission from 
B1800 and its PWN. 

The brighter, compact (w 7" x 4", i.e. 0.14 x 0.08 pc- at 
d = 4 kpc) PWN component is elongated in the southeast- 
northwest direction, roughly perpendicular to the proper mo- 
tion direction (see Fig. 1). The compact component (the iimer 
PWN) is brighter northwest of the pulsar but more extended 
toward southeast. There is also evidence for a more extended, 
fainter emission component (the outer PWN), ^ 12" in size, 
mostly concentrated southward of the compact bright compo- 
nent (Fig. 2). The bottom right panel in Figure 1 shows the 
one-dimensional distribution of counts along the southeast- 
northwest direction. The data points with error bars (con- 
nected by the dashed line) are obtained by integrating counts 
within 10x1 pixels (5" x 0"5) rectangular apertures (shown 
in the top right panel of Fig. 1) moving along the southeast- 
northwest direction with a 0."5 step. The solid line shows the 
count distribution obtained in a similar way but for a point 
source simulated with MARX^. The difference between the 
dashed and solid lines represents the extended emission from 
the PWN. The one-dimensional profile and the PWN images 
suggest a compact (~ 4" - 5" in diameter), roughly synmet- 
ric torus around the pulsar (the blue region in the bottom left 
panel of Fig. 1), with some extension toward southeast. The 
position angle (PA) of the toms symmetry axis is about 48°- 
50° east of north, close to the proper motionPA of 38.°1±6.°3 
(Brisken et al. 2006). 

We also attempted to search for signatures of an SNR 
around the pulsar-PWN complex. A direct visual inspection 
of the ACIS image did not show clear signatures of large-scale 
diffuse emission. We appUed the exposure map correction and 
smoothed the image with various scales, but failed to find sta- 
tistically significant deviations from a uniform brightness dis- 
tribution in these images. To estimate an upper limit on the 
SNR emission, we measured the count rate from the entire S3 
chip (with all identifiable point sources removed). The count 
rate, 0.596 ± 0.004 counts s"' in the 0.5-7 keV band, ex- 
ceeds the nominal S3 background of 0.32 counts s"' (Chandra 
Proposers' Observatory Guide-', v.8, §6.15.2), which could 
be caused by an elevated particle background, diffuse X-ray 
background, or SNR emission. Since we see no trace of an 
SNR, we consider the difference, 0.28 counts s"\ as an upper 
limit on the SNR count rate in the 70 arcmin^ of the chip area, 
which corresponds the average surface brightness limit of 4 
counts ks~' arcmin"-^. 

2.2. Spectral analysis 

2.2.1. PWN spectrum 

The PWN spectra were extracted from two regions shown 
in Figure 1 (top right panel). The smaller (polygon) region 
of 32.7 arcsec^ area encompasses the brighter compact PWN, 
while the larger (circular) region of 176 arcsec^ area includes 

Brisken et al. (2006) position and the position of the X-ray source. No X-ray 
source is seen at or near the radio timing position. This means that the error 
of the radio timing position is strongly underestimated. 

^ MARX (Model of AXAF Response to X-rays) is a suite of pro- 
grams designed to enable the user to simulate the on-orbit perfor- 
mance of the Chandra satellite, including ray-trace simulation. See 
http://space.mit.edu/ASC/MARX/ 

^ See http://asc.harvard.edu/proposer/POG/index.html 
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Fig. 1. — Top left: ACIS-S3 image of B180() and its PWN. Top right: Extraction regions used for the spectral and image analysis of the source components 
(green and white lines; see text). The white arrow shows the direction of the pulsar proper motion (Brisken et al. 2006); its length corresponds to the projected 
distance traveled by the pulsar in 300 years. Bottom left: The smoothed sub-pixel (0.25 of the original ACIS pixel) resolution image obtained by removing the 
pipeline pixel-randomization and subsequently applying the sub-pixel resolution tool, based on analyzing the charge distribution produced by an X-ray event 
(Tsunemi et al. 2001; Mori et al. 2001) and maximum-likelihood deconvolution procedure (Richardson 1972; Lucy 1974). The straight line shows the rotation 
axis of the putative torus (see text). Bottom right: Intensity distribution along the line perpendicular to the rotation axis measured in Of'5 X 5" box moving along 
the line shown in the top right panel. The solid line shows the simulated one-dimensional point-spread function of the Chandra ACIS. 
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Fig. 2.— Adaptively smoothed ACIS-S3 image of B1800 and its PWN. 
The brightness and smoothing scales are chosen to show the fainter, more 
extended emission from the PWN. 



a fainter component, which is more prominent south of the 
pulsar. From these regions we excluded the circular region 
of 1 ."46 radius centered on the brightest pixel to avoid con- 
tamination of the PWN spectrum by the pulsar. The back- 
ground was measured in the 20" < r < 25" annulus centered 



on the source (45 counts in the 706 arcsec area, in the 0.3-8 
keV band). The total numbers of counts extracted from the 
smaller and larger PWN regions are 58 and 117, of which 
96.4% and 90.3 % are expected to come from the source, 
which gives 55.9 ±7.6 and 105.6 ± 11.6 PWN counts in the 
two regions. (The errors here and below are at the 68% confi- 
dence level.) The observed PWN fluxes are Fp^n = (5.5 ± 0.6) 
and (2.7 ± 0.4) x lO"'** ergs s"' cm"^ for the larger and smaller 
extraction regions, respectively, in the 1-8 keV band. 

For each of the two spectra, we group the counts into six 
spectral bins, with comparable numbers of counts per bin. 
First, we fit the PWN spectra for each of the two regions with 
the absorbed PL model, allowing the hydrogen column den- 
sity, tin, to vary. These fits result in spectral slopes Fpwn ~ 
1.5-1.9 and «h,22 = n^/lO^^-cmr^ « 1.2-2.0 (see Fig. 3). We 
see no indication for spectral softening (expected due to syn- 
chrotron burn-off in PWNe with high radiation efficiency) in 
the spectrum extracted from the larger region compared to the 
spectrum of the bright inner PWN. To test this, we also fit 
the spectra with the PL model in which the hydrogen column 
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Fig. 3. — Confidence contours (68% and 90%) in tlie n}\-V plane for the 
PL fit to the PWN spectrum. The error bars show the photon indexes for the 
fixed «H.22 = 0.72. Blaclc and green contours and error bars correspond to the 
larger and smaller PWN extraction regions, respectively (see text and Fig. 2). 

density is fixed at nH,22 = 1-38 (the best-fit value for the larger 
region) and found that the virtually the same photon indices: 
Tpwn = 1 .58 ± 0.25 and 1 .55 ± 0.29, for the larger and smaller 
region, respectively. The apparent lack of spectral softening 
could be attributed to the small number of counts collected 
from the fainter outer PWN component. 

The best-fit value of n-^xi ~ 1-4 estimated from the B1800 
PWN fit approximately coincides with the total Galactic HI 
column density in this direction (Dickey & Lockman 1990). 
It does not contradict to the adopted distance of «4 kpc to 
B1800 because the hh deduced from an X-ray spectrum un- 
der the assumption of standard element abundances generally 
exceeds the hhi measured from 21 cm observations by a factor 
of 1 .5-3 (e.g., Baumgartner & Mushotzky 2005). On the other 
hand, given the B1800's dispersion measure, DM = 234 cm"^ 
pc (i.e., the electron column density He = 7.22 x 10^° cm"^), 
the nH,22 value of 1.4 corresponds to the ISM ionization de- 
gree Hc/nH ~ 5%, below the usually assumed value of 10%, 
which corresponds to nH,22 = 0.72. If the hydrogen column 
density is fixed at nH.22 = 0.72, the PL fit results in substan- 
tially harder spectra, Fpwn ~ 1 .0 for both extraction regions 
(see Table 1 and Figs. 3 and 4). 

The unabsorbed flux and luminosity are less sensitive to 
the interstellar absorption than the spectral parameters. The 
isotropic PWN luminosity is Lpwn = 47rii^Fp"™'" « 1-2 and 

0.6-1 X lO^^dl ergs s"' in the 0.5-8 keV band, for the whole 
and inner PWN, respectively (see Table 1 and Fig. 4). The for- 
mer value is a fraction of « (0.5-1) x lO'^^ii^ of the pulsar's 
spin-down power £. 

2.2.2. Pulsar spectrum 

To minimize contamination by the PWN, the pulsar spec- 
trum was extracted from a small circular aperture with the 
radius of 1.5 ACIS pixels (~ 0."74, 85% encircled energy ra- 
dius). The total number of counts within the source region is 
5 1 , of which less than 1 5 % is expected to come from the PWN 
(based on the PSF simulation). To account for the PWN con- 
tribution to the pulsar spectrum, the background spectrum was 
extracted from the 8 arcsec^ region encompassing the bright- 
est part of the inner PWN (with the pulsar being excluded). 
The pulsar's absorbed flux is Fpsi- = (1.4 ±0.2) x 10"'"^ ergs 
cm"^ s"' in the 1-6 keV band (aperture corrected and back- 



FIG. 4.— Confidence contours (68%) and 90%) for the PL fit to the PWN 
spectra with nH,22 = 1-38 (black) and 0.72 (red). The upper contours corre- 
spond to the larger extraction region. The PL normalization is in units of 10"^ 
photons cm"^ s"' keV"' at 1 keV. The dashed curves are the loci of constant 
unabsorbed flux in the 0.5-8 keV band; the flux values near the curves are in 
units of 10"'* ergs cm"^ s"' . 
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Fig. 5. — Confidence contours (68% and 90%) for PL spectral parameters 
for various fits to the pulsar spectrum: PL fit with 22 = 0.72 (blue). PL+BB 
fit with nH.22 = 1-38 (black), and PL+BB fit with hh.22 = 0.72 (red). The PL 
normalization is in units of 10"* photons cm 2 s"' keV"' at 1 keV. The dashed 
curves are the lines of constant unabsorbed flux in the 0.5-8 keV band (the 



flux values are in units of 10 ergs cm ^ 1 



ground subtracted). To obtain constrained fits with the small 
number of counts available, we are forced to freeze the hy- 
drogen column density. First, we fix it at «h,22 = 1-38, ob- 
tained above from the PL fit to the PWN spectrum. With 
this «H, the single-component PL and black-body (BB) mod- 
els fail to fit the pulsar spectrum. A two-component BB-i-PL 
model provides a good fit, but the fitting parameters are poorly 
constrained because of the small number of photons detected 
(Fig. 5). The slope of the PL component is Fpsr = 1.4 ±0.6, 
and its unabsorbed luminosity is L™"* = (4.4 ± 1 . 1) x 10^'^^| 

ergs s-> « 2 X lO'^Edj, in the 0.5-8 keV band. The tem- 
perature and the projected area of the BB component are 
strongly correlated (see Fig. 6), which results in very large 
uncertainties for these parameters. The best-fit temperature is 
r w 1 .6 MK, while the projected emitting area, A^2x 10^c/| 
m^, is smaller than that of the NS surface (nR^ ~ 3 x 10^ 
m^), but larger than the conventional polar cap area Ape = 
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Fig. 6. — Confidence contours (68% and 90%) for the BB component of the 
BB+PL fit to the pulsar's spectrum, for )7h.22 = 1-38 (black), and jih,22 = 0.72 
(red). The BB normalization (vertical axis) is the projected emitting area in 
units of m^, assuming the distance of 4 kpc. The lines of constant bolometric 
flux (in units of lO"'* ergs cm"^ s"') are plotted as dashed lines. 



TABLE 1 
PL FITS TO THE PWN SPECTRUM 



Model 


nH.22" 




r 


(C or x^)'/dof 




Entire PWN 
Entire PWN 
Inner PWN 
Inner PWN 


1.38 
0.72 
1.38 
0.72 


13.8+3-5 
6.0:11 

9.7:3.7 

3 3+' 2 


1 co+0.25 
^•^''-0.25 
1 1+0-23 

1 70+0.29 
-^-0.32 

1 04+0-29 
^•"^-0.26 


0.93/3 
2.16/4 
1.34/3 
1.34/4 


1.65!!!;15 

1 35+0-" 

96+'"6 
-0.13 

70+" " 
-0.10 



"The values of hydrogen column density are frozen in the fits. 

''Spectral flux in units of 10"* photons cm"2 g-l keV"' at 1 keV. 

■^We use the C statistic (Cash 1979) for the inner PWN and statistic for 
entire PWN, which has a larger number of counts. 

''Unabsorbed isotropic luminosity in the 0.5-8 keV band in units of lO'^rfl 
ergs s"' . 

l-K^F? jcP « 5 X 10^ m^. The corresponding bolometric lu- 
minosity, LJ^°' = ^AoT^ - 3 X lO^^t/l ergs s"'. 

If the hydrogen column is fixed at the value estimated from 
the dispersion measure at 10% ISM ionization, «h,22 = 0.72, 
then a single PL model fits the spectrum and gives Fpsi = 
2.0!;; j and the unabsorbed luminosity L™"* w (4.9 ±0.8) x 
lO^'c/4 ergs s"', in the 0.5-8 keV band. Once again, the BB 
model does not fit the spectrum while the two-component 
PLh-BB model provides an acceptable fit. In comparison with 
the BBh-PL fit with the larger nH.22 = 1-38, the BB compo- 
nent shows a higher temperature ij « 2.6 MK) and a smaller 
(but even more uncertain) emitting area (^ 2 x 10^t/| m^) 
while the PL component has Fpsi « 1.2. The fitting pa- 
rameters of the BB component are even less constrained 
than those at nH,22 = 1-38 (see Fig. 6). For the non-thermal 
PL component, the confidence contours for fitting parame- 
ters strongly overlap with those calculated for nH.22 = 1-38 
(see Fig. 5). The luminosity of the non-thermal component, 
^nonth ^ (3.7 ±0.9) X 10^'t/| ergs s"' in 0.5-8 keV, is close 
to that in the n//,22 = 1 -38 case, while the bolometric luminos- 
ity of the thermal component, Lp°' ~ 2 x 10-"t/| ergs s"', is a 
factor of 15 lower 

3. DISCUSSION. 

The B1800 PWN-npulsar flux measured by Chandra is con- 
sistent with the ROSAT PSPC count rate reported by Finley & 
Ogelman (1994). However, ROSAT was unable to resolve the 



PWN emission from that of the pulsar. The superior angular 
resolution of Chandra has allowed us to resolve the shape of 
the PWN, disentangle the point source and the extended emis- 
sion components, and study their properties separately. Below 
we discuss these properties in more detail and compare the re- 
sults on B 1800 with those on other Vela-like pulsars. 

3.1. Energetics and spectra ofBlSOO and other Vela-like 
pulsars and PWNe 

Since the pulsar wind energetics is supplied by the loss of 
the pulsar spin energy, it is natural to expect that the PWN 
luminosity is correlated with the pulsar spin-down power, E. 
Such a correlation was first noticed in X-ray PWN observa- 
tions with the Einstein observatory (Seward & Wang 1988). 
On the other hand, as the wind properties can depend on, e.g., 
pulsar's magnetic field and the angle between its magnetic 
and spin axes, and the PWN properties depend on pulsar's ve- 
locity and pressure of ambient medium, we should not expect 
that the X-ray PWN efficiency, rjp^n = Lpvm/E, is the same 
for all PWNe. To examine the Lp^„-E correlation for similar 
PSR/PWN systems (t = 10-30 kyrs, E = 10-^''-10" ergs s-^)\ 
we have estimated the X-ray PWN luminosities in the 0.5- 
8 keV band for 9 such objects (in addition to B1800) using 
archival Chandra data (see Table 2). Figure 7 demonstrates 
that although there is a positive correlation between Lp„n and 
E, it is far from linear and shows a large scatter In particu- 
lar, the PWN efficiency spans over 2.5 orders of magnitude, 

10-4-5 < 

^pwn ^ 10 ^ in the relatively narrow E range, (2— 
7) X lO^*" ergs s"' (the B 1800 PWN luminosity and efficiency 
are among the lowest in this sample). Such a scatter can 
hardly be explained by poorly known distances to some of the 
objects or by the bias caused by a possible underestimation of 
Lpvm for distant PWNe, where we can miss a component with 
lower surface brightness. Interestingly, the two bow-shock 
PWNe with prominent tails, the Mouse and the Duck, show 
efficiencies higher than most of the other PWNe in the sample, 
perhaps because the wind is channeled into a smaller volume 
produces a brighter X-ray image. Another very bright and ef- 
ficient PWN is that around PSR J181 1-1925, which is located 
in the young SNR Gl 1.2-0.3, possibly associated with the 
historical SN of 386 AD (Kaspi et al. 2001b). If this associa- 
tion is correct, the true pulsar age is a factor of 14 smaller than 
the spin-down age, so that PSR J181 1-1925 is much younger 
than the other pulsars in the sample, which may explain why 
its PWN is so different. Excluding these three PWNe from the 
sample substantially narrows the ranges of PWN luminosities 
and efficiencies: 

Vn-^lO^'-lO^-^ergss-i, 7?pw„ ^ lO'^'^ - 10"-^ ^ (1) 

It is also interesting to compare the PWN X-ray luminosi- 
ties with the magnetospheric (PL component) pulsar lumi- 
nosities. We see from Figure 7 that the former exceeds the 
latter for 9 out of the 10 objects in our sample (the only ex- 
ception is the Duck PWN and its pulsar B 1757-24). More- 
over, we see a clear correlation between the two luminosities, 
which can be described as logLp„n = 0.89 logL™"* -1-4.2 or, 
more crudely, 

r ^rnonth 

in the 0.5-8 keV band. Overall, the correlation is surpris- 
ingly strong, given the fact that Lp°"* was estimated assum- 

4 Pulsars with such t and E are traditionally dubbed Vela-like pulsars. 
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Fig. 7. — PWN luminosity versus pulsar non-thermal luminosity (top) 
and spin-down power (bottom) for 10 Vela-like pulsars observed with Chan- 
dra. The luminosities are estimated for the 0.5-8 keV band. The solid and 
dashed lines in the top panel correspond to logLpwn - 0.89 log Lp™"" +4.2 and 

Lpwn = 5Z,p"*, respectively. The dash-dot Unes in the bottom panel are lines 
of constant PWN efficiency, rypwn- The error bars include the statistical un- 
certainties (see Table 2) and the nominal 30% distance uncertainties, except 
for the Vela pulsar whose parallax has been measured (Dodson et al. 2003). 

ing isotropic pulsar emission while the actual pulsar lumi- 
nosity can be higher or lower than this estimate. The main 
outliers from this dependence are the Vela, whose pulsar 
is unusually underluminous (Pavlov et al. 2001a), and the 
Duck/B 1757-24, for which separation of the pulsar and PWN 
emission is rather unreUable even with the Chandra resolution 
(Kaspi et al. 2001a). On the other hand, B1800, for which 
Lpwn/ip™* ~ 4, is a typical PSR/PWN pair in this regard. A 
natural reason for such a correlation is that both the magne- 
tospheric and PWN radiation are powered by the spin-down 
energy loss (i.e. both Lp„n and Lp°"* are correlated with E and 
hence with each other). However, the Lp^n-i^ps"* correlation 

shows a smaller scatter than the L^m-E correlation, despite 
the quite different nature of the magnetospheric and PWN X- 
ray emission. This may suggest that the L^vm-E correlation 
is significantly distorted by the distance errors (in contrast to 
the Lp^4,n-Lp°"* correlation that does not depend on distance) 
or, more likely, that there are some other factors (e.g., the an- 
gle between the spin and magnetic axes) that similarly affect 



the magnetospheric and PWN emission. 

The spectral slope of the B1800 PWN is somewhat uncer- 
tain because of the small number of counts and the depen- 
dence of the fitting parameter Tp^n on the poorly known mh 
(see Figs. 3 and 4). For a plausible nH,22 = 1-4, the pho- 
ton index, Fpwn — 1.6 ±0.3 for the entire PWN, is typical 

for the whole sample of PWNe observed in X-rays^, and 

'-ill 

it is substantially larger than Fpwn = 2.36-2.1£'3g' = 0.94 
(where £36 = £/10^^ergss"') predicted by the Tp^-E corre- 
lation found by Gotthelf (2003) for a sample of more ener- 
getic pulsars. On the other hand, a lower assumed nH,22 = 0.7, 
which cannot be ruled out based on the data available, results 
in an unusually hard PWN spectrum, Fp^n ~ 1.0 ±0.2. For the 
broad range of ne considered, Fpwn remains similar to those 
of the other PWNe in Table 2 except for the Mouse and the 
Duck, whose softer spectra (Fp^n ~ 2.0-2.5) can possibly be 
explained by a stronger effect of synchrotron cooUng in the 
long PWN tails. 

The spectrum of the B1800 pulsar is even more uncertain 
than that of its PWN. As we have mentioned in §2.2.2, the 
one-component PL model does not fit the pulsar spectrum for 
«H,22 = 1 -4, while the PL component of the PL-hBB fit gives 
Fpsr ~ 1.4 ±0.6, similar Fp™ of B1800 and to the spectral 
slopes of many other pulsars. This slope is much softer than 

Fpsr = l.l-1.9Ell'^ = 0.1 predicted by the Gotthelf 's (2003) 
correlation. Lowering the assumed «h by a factor of 2 makes 
a single-component PL fit acceptable, resulting in a rather soft 
spectrum, Fpsr — 2.0 ± 0.4. The PLh-BB fit at this smaller tin 
yields a harder PL component, with Fpsr « 1 .2 (see Fig. 5), 
which is close to the Fp™ obtained for this rtn, but still con- 
siderably softer the Gotthelf 's (2003) prediction. Overall, we 
can only state that Fpsr for B1800 is within the range of val- 
ues observed for the whole sample of radio pulsars detected 
in X-rays, including the Vela-hke pulsars. 

The possible thermal component of the PLh-BB fit to the 
B1800 pulsar spectrum is very poorly constrained, not only 
because of the poor statistics but also because the soft thermal 
radiation is strongly absorbed by the ISM. The BB tempera- 
ture, T ~ 1-3 MK, and bolometric luminosity, L^l 10^^- 

10^^ ergs s"\ are similar to those found from the PLh-BB fits 
for other Vela-like pulsars. Such temperatures are somewhat 
higher than the bulk NS surface temperatures predicted by the 
standard NS cooling models (e.g., Yakovlev & Pethick 2004), 
and the corresponding emitting areas are smaller than the NS 
surface area. However, the actual spectrum of the NS ther- 
mal radiation can differ substantially from the BB model. In 
particular, fitting the spectra with the hydrogen atmosphere 
models (Pavlov et al. 1995) yields lower effective tempera- 
tures and larger emitting areas, with not so strongly different 
bolometric luminosities (see Pavlov et al. 2001a for the spe- 
cific example of the Vela pulsar). Unfortunately, the quality 
of the data do not warrant fits with more complicated atmo- 
sphere models, and even the estimate for the bolometric lumi- 
nosity is too uncertain to make a useful comparison with the 
NS cooling models. 

3.2. Inferences from the PWN morphology 



We should caution that the photon index generally grows with increasing 
distance from the pulsar due the synchrotron cooling. Therefore, the Tpwii 
value depends on the extraction area chosen, and it may be systematically 
smaller for more distant and fainter PWNe, where the outer regions with 
softer emission are too faint to include them in the spectral analysis. This 
observational bias must be taken into account when comparing PWN spectra. 
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Fig. 8.— Chandra ACIS image of the Vela PWN in tlie 1-8 keV band as it 
would be seen at the distance of 4 kpc in a 30 ks exposure. The resolution is 
degraded by a factor 1 3 to simulate the Vela PWN appearance at the distance 
of 4 kpc. The aiTow shows the direction of the Vela pulsar proper motion. 



The distribution of X-ray counts in the immediate vicinity 
of the B1800 pulsar shows that the X-ray PWN is elongated 
in the direction perpendicular to the pulsar's proper motion. 
By analogy with young, bright PWNe (such as the Crab and 
Vela), the observed PWN shape can be interpreted as a torus 
with a symmetry axis along the southwest-northeast direction 
(see Fig. 1). The torus appears to be seen nearly edge-on, 
implying that the pulsar's spin axis is close to the plane of the 
sky. We find no evidence for jets along the symmetry axis. 
This, however, should not be surprising because such jets can 
be significantly fainter than the torus (for instance, the Vela 
PWN outer jets would not be seen in a 30 ks exposure if the 
Vela PWN were located at a 4 kpc distance while the inner 
jets would not be resolved; see Fig. 8). 

In the B1800 PWN, the torus symmetry axis, which pre- 
sumably coincides with the pulsar spin axis, is approximately 
aligned (within « 10°) with the direction of the pulsar's 
proper motion. This strengthens the evidence for the align- 
ment between the pulsar spin and velocity vectors, based 
on a very limited sample of pulsars so far (e.g., Ng & 
Romani 2004; Johnston et al. 2005), which constrains the 
physics of SN explosions. It also implies that the total (three- 
dimensional) speed of the pulsar is close to the measured 
transverse velocity in the plane of the sky, vj^ = (365 ± 30)ii4 
km s"' (Brisken et al. 2006). 

Toroidal (or arc-like) X-ray structures are often found 
around young pulsars, the Vela and the Crab being the most 
famous examples. Out of the ten Vela-like pulsars observed 
with Chandra (see Table 2 and Fig. 9), three (Vela, B1706- 
44, and J202 1-1-3651) clearly show arc-like (or toroidal) struc- 
tures and jets in X-rays. On the other hand, the PWNe around 
J1747-2958 (the Mouse) and B1757-24 (the Duck) exhibit 
different X-ray morphologies, with prominent tails behind 
the moving pulsars, indicative of bow-shock PWNe. For the 
remaining four pulsars (B1823-13, J1016-5857, B1046-58, 
and J1811-1925), the existing X-ray data do not allow one to 
firmly claim the existence of either a tail or a torus based on 
X-ray data alone. 

Despite the relatively high speed of B 1800, the PWN shows 
no evidence of bow-shock morphology. This may look sur- 
prising since, for instance, the Mouse PWN, whose pulsar has 
almost the same age and spin-down power as B1800, clearly 
shows such morphology, although the speed of its pulsar, in- 



ferred from the bow-shock modeling, is only a factor of 1.6 
higher (Gaensler et al. 2004)^. Furthermore, for the Duck 
PWN, recent radio-interferometric measurements put an up- 
per limit of 340((i/5kpc) km s"' on the transverse pulsar 
speed (Blazek et al. 2006). Nevertheless, the Duck shows a 
distinct trail in both radio and X-rays (Frail & Kulkarni 1991; 
Kaspi et al. 2001a), which is a signature of a ram-pressure 
confined, supersonically moving PWN. Thus, the different ap- 
pearance of the B 1 800 PWN as compared to the Duck and the 
Mouse should be attributed to differences in the properties of 
ambient medium rather than to different pulsar velocities. 

Assuming that the lack of bow-shock morphology in B 1 800 
implies a subsonic pulsar velocity, we obtain a lower limit 
on the temperature of the ambient medium, T > lOdj MK. 
Such a high temperature can only be expected inside a young 
SNR (e.g., temperatures of up to 14 MK have been measured 
in the vicinity of the Vela pulsar, although typical tempera- 
tures in the Vela SNR are 1-3 MK; Lu & Aschenbach 2000). 
This means that the radius of the putative host SNR is sub- 
stantially larger than v^r « 6ii4 pc (/? > 5' in angular scale; 
for comparison, the radius of the Vela SNR, ^ 4° at 300 pc, 
would be ^ 18' at d = 4 kpc). However, our analysis of the 
ACIS image showed no evidence of an SNR (see §3.1), nei- 
ther the putative host SNR of B1800 nor the G8. 7-0.1 rem- 
nant onto which B1800 is projected. If the X-ray spectrum 
and luminosity of the undected SNR are similar to those of 
the Vela SNR, the nondetection is not surprising because of 
the much larger nu for B 1 800. For instance, using the PIMMS 
tool^, we estimate that the average unabsorbed intensity from 
the Vela SNR, I ^ 1.5 x lO"'"* ergs cm"^ s"' arcmin"^ in the 

0. 1-2.4 keV band (Lu & Aschenbach 2000), corresponds to 
the on-axis ACIS-S3 surface brightness of < 0.3 counts ks"' 
arcmin"^ in the 0.5-7 keV band, for the Raymond-Smith ther- 
mal plasma emission models with T <3 MK and nu.22 = 1 -4, 

1. e. at least a factor of 13 below the observed upper limit. 
Even the hot thermal component of the Vela SNR, such as ob- 
served in the vicinity of the Vela pulsar, could hardly be de- 
tected in our observation if Vela were at the location of B 1 800. 
For instance, / ^ 3 x 10"'^ ergs cm"^ s"' arcmin"^ in the 0.9- 
2.4 band, as estimated from Figures 15 and 9 of Lu & As- 
chenbach (2000), corresponds to the ACIS-S3 surface bright- 
ness of only 1.2 counts ks"' arcmin"^, for T = 10 MK and 
nH.22 = 1-4. On the other hand, Finley & Ogelman (1994) did 
detect extended X-ray emission ^ 30' northeast of the pulsar 
(this region was out of the Chandra ACIS field of view) and 
interpreted it as due to a shock in the G8. 7-0.1 SNR, reflected 
from a nearby molecular cloud. Fitting the spectrum of that 
emission with the Raymond-Smith model, they found T = 4-8 
MK and 11^ 22 = 1.2-1.4. Based on the currently availble data, 
we cannot conclude whether this emission is associated with 
the putative host SNR of B1800 or with G8.7-0.1. In the lat- 
ter case, the G8. 7-0.1 SNR should be at a distance similar to 
that of the host SNR of B1800 (e.g., '-1 kpc foreground) as 
follows from the similar «h values. 

* We should note that the proper motion of PSR J1747-2958 has not been 
measured. The Mach number Ad = v/cj = 60 was obtained by Gaensler at 
al. (2004) assuming an ambient pressure Pamb = 3.3 X 10"'^ ergs cm~^, and 
the pulsar speed v = 600 km s"' was estimated assuming that the ambient 
medium is a warm ISM with the sound speed c, « 100 km s"' . Both M and 
I' would be smaller if the ambient pressure is higher than assumed, as one can 
expect for such a young pulsar (r = 26 kyr), which is possibly at outskirts of 
its SNR. For instance, X = 25(rf/5kpcr' for 

Pdmh — 2 X 10 " ergs cm ^ , and 
V = 355«"'^^(rf/5kpc)"' km s"', close to the transverse velocity of B1800. 
' See http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html. 
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Fig. 9. — Chandra ACIS images of PWNe around Vela-like pulsars (see Table 2). 
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The shape and size of the B1800 PWN can be used to ob- 
tain additional constraints on the properties of the ambient 
medium, such as pressure, density, and magnetic field. For 
instance, the ambient pressure confining the PWN should be 
a sizeable fraction of the pulsar wind pressure at the termina- 
tion shock, pamb = fpPs, where fp<\, 

Ps = EMA^crl)-\ (3) 

Vs is the termination shock radius in the equatorial plane, and 
the factor /q takes into account anisotropy of the pulsar wind 
(for instance, /o ~ 3/2 for the often assumed dependence of 
the wind flux, oc sin^^, on rotational colatitude 6; e.g., 
Bogovalov & Khangoulyan 2002). In well-resolved bright 
PWNe, such as the Crab and Vela, we see the "inner ring", 
commonly associated with the pulsar wind termination shock 
in the equatorial outflow. We cannot resolve such a ring in the 
B1800 image, but, by analogy with the bright, well-resolved 
PWNe, we can assume its radius to be about one half of the ra- 
dial extent of the X-ray bright torus in equatorial plane. Based 
on the surface brightness profile shown in Figure 1 (bottom 
right panel), we take rj « 1"5 « O.Q2>d4 pc as a reasonable 
estimate, which is close to the values measured in other 
Vela-like X-ray PWNe with clearly visible toroidal structures 
(including Vela, B 1706^4, and J2021H-3651; see Table 2). 
Then, the ambient pressure can be estimated as 

Pamb ~ 7.2 X IQ-^'^fpfndt ergs cm'^. (4) 

Similar to many other PWNe around young pulsars (see Fig. 
10), the estimated ambient pressure is much higher than the 
typical unperturbed ISM pressure, ^ 10"'^ ergs cm~^, in 
agreement with our inference that B 1 800 has not left its SNR. 

Using a similar approach, we estimated the pulsar wind 
pressure at the termination shock for 13 other young PWNe, 
including the Mouse and the Duck for which is the es- 
timated separation between the pulsar and the bow-shock 
head, and plotted ps versus pulsar spin-down age in Fig- 
ure 10, assuming /n = 1- Even with allowance for large 
uncertainties, the pressure does not follow the dependence 
Pamb c< T~^l^ expected for SNRs in Sedov stage, nor any 
other smooth dependence. Particularly surprising is the low 
pressure, ;7j ~ 4 x 10~'^/n(<i/4.4kpc)~'^ ergs cm"^, for the 
very young (r ~ 1.6 kyr) PSR B1509-58, which suggests a 
very unusual SNR for this object, in agreement with X-ray 
observations (Gaensler et al. 2002). It also hints that PSR 
B 1509-58 moves supersonically (which would require a ve- 
locity above 50(t//4.4kpc)~^n"'/^ km s~\ where n = p/niH), 
so that the arc(s) northwest of the pulsar and the bright, long 
"jet" southeast of the pulsar are elements of a bow-shock 
PWN rather than a "torus-jet" PWN. We should also note that 
the spin-down age of a pulsar can differ substantially from 
its true age, which can contribute to the lack of a clear p^-r 
correlation in Figure 10. For instance, r = 620 kyr for PSR 
J0538H-2817 is a factor of 6 larger than the estimated age of 
the G180.1-1.7 SNR, presumably associated with this pulsar 
(e.g., Romani & Ng 2003). 

Interestingly, the pulsar wind pressure at the head of bow- 
shock PWNe is similar to that at the termination shock radius 
of torus-like PWNe, despite the different confining mecha- 
nisms. The reason for this coincidence is that, for typical pul- 
sar speeds, the ram pressure caused by the pulsar motion in 
file ISM, /7ran, = pv^ = 1 .5 X 10-''n(v/300kms"')2 ergs cm'^ 
is of the same order of magnitude as the ambient pressure in- 
side an SNR confining a torus-like PWN. 



Using the above estimate for the ambient pressure confining 
B 1 800 and assuming that the ram pressure due to the pulsar's 
motion, 

Pram » 2.2 X lO'^ndj ergs cm"^ , (5) 

is lower than pamb (as expected for the subsonic motion), we 
obtain an upper limit on the density of the ambient medium, 
n < 03fpfndf cm-3, or p < 5 x IQ-^^fpfnd^ g cm'^. 

The measured termination shock radius can also be used to 
estimate the magnetic field inside the X-ray PWN. According 
to Kennel & Coroniti (1984), the particle pressure immedi- 
ately downstream of a strong termination shock perpendicular 
to the wind (i.e., close to the equatorial plane in our case) is 
2ps/3, so that the internal energy is 2ps. If cb is the fraction 
of the internal energy in the magnetic field, the field can be 
estimated as B = (167reB/>,)'/^, or 

B=(l!^y'^190ey^/^/VMG. (6) 

The fraction eg varies along the post-shock flow, and it de- 
pends on the magnetization and angular distribution of the 
pre-shock pulsar wind, as well as on the properties of the 
termination shock. For instance, in the Kennel & Coroniti 
(1984) model of spherical shock in an isotropic wind with a 
small pre-shock magnetization parameter a, the magnetic en- 
ergy fraction is ~ 9c7/4 immediately downstream of the 
shock, and it grows with increasing radius up to ~ 1 at 
r w r,(3<T)-'/2. In 

a more realistic equatorial outflow geome- 
try, Komissarov & Lyubarsky (2004) found a highly nonuni- 
form distribution of eg in the postshock flow, with maximum 
values close to unity. Given the uncertainty of the flow pat- 
tern in B1800, we can only crudely estimate ~ 100 /xG for 
a typical magnetic field in its PWN, similar to the Vela and 
other Vela-like PWNe (e.g., Gonzalez et al. 2006). 

A crude estimate of the synchrotron cooling time for elec- 
trons emitting at an energy E (in keV) is Tsyn ~ 36BI^^^£'"'/^ 
yrs, where B_4 = B/lOO/iG. During the time Tsyn, the pulsar 

travels the projected distance, v_LTsyn « 0.0l3d4E~^/^B^^^ pc, 

corresponding to the angular distance of 0'.'7E~^^^B^^^, con- 
siderably smaller than the size of the compact X-ray PWN. 
For a typical shocked wind flow velocity of ~ 0.3c, the elec- 
trons travel a distance of 33E^^^^B'^^^ pc during the syn- 
chrotron cooling time, corresponding to an angular distance 

of 2'.Sd4E~^^^B'^^'^, which substantially exceeds the observed 
extent of the X-ray PWN. We caution, however, that the flow 
streamlines are curved for an anisotropic pulsar wind (see, 
e.g., Komissarov & Lyubarsky 2004), so that the actual PWN 
size can be smaller than this distance. In addition, the syn- 
chrotron surface brightness decreases with the distance from 
the pulsar (see e.g., the Vela PWN image in Fig. 9) rendering 
fainter emission undetectable during the relatively short ACIS 
exposure. This could also explain the apparent lack of spec- 
tral softening (§2.2.1) due to the synchrotron burn-off, which 
should become more significant further away from the pulsar. 
A deeper exposure can reveal the expected faint extension of 
the X-ray torus beyond the 5" radius if the equatorial flow 
remains sufficiently well collimated. 

The compact inner PWN of B1800, which we interpret as 
a torus seen edge-on, shows some asymmetry with respect 
to the spin axis: the projected torus is more elongated to- 
wards the southeast, and its southeast edge is more diffuse 
and fainter than the northwest edge (see Fig. 1). A possi- 
ble explanation for this asymmetry is a pressure gradient in 
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Log Age, yrs 

Fig. 10. — Pulsar wind pressure at termination shock versus pulsar spin- 
down age for a sample of 14 pulsars, assuming an isotropic pulsar wind. The 
diamonds mark PWNe showing torus-like structure (B1800 is shown by the 
large filled diamond), while the two PWNe with a clear bow-shock struc- 
ture are shown by squares. The pulsars B1800-21, B1757-24 (the Duck), 
J2021+3651, and B 1046-58 are denoted by the numbers 1, 2, 3 and 4, respec- 
tively. To calculate the shock radius r,, we used the distance d in the sixth 
column of Table 2. The straight line, arbitrarily drawn through the B1800 
point, shows the dependence p tx r~*/' expected for an SNR in Sedov stage. 

the ambient SNR matter. As r, oc /J^nib' pressure gradient 
leads to an azimuthal dependence of the termination shock ra- 
dius in the equatorial plane, such that the circular ring turns 
into an ellipse-like structure. Since the surface of such an 
asymmetric shock is not perpendicular to the unshocked pul- 
sar wind even in the equatorial plane, the post-shock flow is 
deflected toward the most distant (southeast) part of the shock, 
where the magnetic field and the synchrotron emissivity are 
lower. To produce the observed asymmetry, the pressure 
difference should be a sizable fraction of the mean ambient 
pressure, Ap ^ a few x 10"'" ergs cm"-'. Such a difference 
should cause a local wind in the ambient medium, possibly 
on a parsec scale, with a typical velocity Vamb ~ i^p/p)^^^ ~ 
240(A/?/3 X 10""'ergscm"3)i/2(p/5 x 10"25g cm"^)"'/^ km 
s"', comparable to the pulsar's velocity. This wind, blow- 
ing, e.g., from the west with a speed comparable to that of the 
pulsar, could also explain the asymmetry of the fainter PWN 
emission, extended toward southeast- south (see Fig. 2). Sim- 
ilar asynmietry is seen in the brighter nearby Vela PWN (see 



Fig. 9), which has also been attributed to a local SNR wind 
(Pavlov et al. 2003). Such pressure gradients and accompany- 
ing winds could also be responsible for the asymmetric shape 
of other PWNe, such as B1823-13 (see Fig. 9). 



4. CONCLUSIONS 

We have detected flie X-ray emission from PSR 61800-21 
and its synchrotron nebula, elongated perpendicular to the di- 
rection of the pulsar's proper motion. The shape of the ob- 
served PWN suggests that, similar to many other young pul- 
sars, B1800 is surrounded by a compact X-ray torus around 
an equatorial termination shock in the pulsar wind. The sym- 
metry axis of the torus (i.e. the pulsar's spin axis) is ap- 
proximately aligned with the direction of the proper motion. 
The compact nebula is embedded into a fainter, more ex- 
tended emission, which shows substantial asymmetry, pos- 
sibly caused by a pressure gradient in the ambient medium. 
Despite a rather high velocity of B 1 800, we find no evidence 
of a bow-shock X-ray nebula. This implies that the pulsar is 
still moving in a hot SNR medium, but we do not see the SNR 
in the Chandra image. In terms of its X-ray efficiency and 
spectrum, the B1800 PWN is very similar to PWNe around 
other pulsars of similar ages and spin-down powers, subson- 
ically moving within their SNRs. Bow-shock PWNe associ- 
ated with young, supersonically moving pulsars show higher 
X-ray efficiencies and softer spectra. 

The spectrum of the B1800 pulsar can be described by a 
two-component BB-hPL model. For reasonable hh values, the 
parameters of this model and the corresponding component 
luminosities resemble those of other Vela-like pulsars. 

A deeper Chandra observation would allow one to inves- 
tigate the morphology of a fainter outer PWN and perhaps 
reveal the jets along the torus symmetry axis. To look for 
the putative host SNR and investigate the PWN emission on a 
larger spatial scale, including its possible connection with the 
TeV source HESS J 1804-2 16, this field should be observed 
with XMM-Newton. 
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TABLE 2 

Properties of Vela-like pulsars observed with Chandra 
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1.5 ±0.3 
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Note. — To ensure uniformity of analysis, we use the results of our own analysis of the Chandra data. These results may differ from the published ones 
because of different extraction regions used, different nn assumed, etc. The uncertainties of the liuninosities and photon indices include the statistical uncertainties 
from spectral fitting but not the systematic uncertainties which can be substantial in some cases and are difficult to calculate. 
"Spin-down power in units of 10^* ergs s"' . 

''Dispersion measure distance from the Cordes & Lazio (2002) model. 

"^Our best guess for a plausible distance to the pulsar, used to scale all the distance-dependent parameters in this table. For three pulsars, J2021+3651, J1016- 
5857, and J1747-2958, the adopted distances differ substantially from dou- In the two cases with the largest DM distances our distance estimates are based on 
the measured nn and the Galactic HI column. For the detailed discussion of the distance to J1747-2958, see Gaensler et al. (2004). 

The hydrogen column density is obtained from the spectral fits to the PWN spectra, except for J1016-5857 and B 1046-58 for which the small number of PWN 
counts precludes a reliable mh measurement. In these two cases the mh values (shown in square brackets) are estimated from the pulsar's DM (assuming 10% 
ISM ionization). For the Vela PWN we quote the luminosity of the inner compact PWN restricted to the "arcs" region (i.e. within 25" from the pulsar). 

'Nonthermal luminosity of the pulsar in the 0.5-8 keV band, in units of lO'- ergs s"' . In the cases when the spectnun is fitted with the BB-hPL model, it is the 
luminosity of the PL component only. 

^Thermal bolometric luminosity of the BB component from the PL+BB fit in units of 10^^ ergs s-'. 

''The largest linear extent of the observed X-ray emission from a PWN. 

j Estimated stand-off distance from the pulsar to the termination shock. 

'^Estimated pressure at the termination shock assuming /n = 1 (see text). 

'References to papers where the corresponding Chandra data have been analyzed. - (1) Pavlov et al. 2001ab; (2) Roberts et al. 2003; (3) Romani et al. 2005; 
(4) Hessels et al. 2004; (5); Teter et al. 2003; (6) Camilo et al. 2004; (7) Kaspi et al. 2001a; (8) Gaensler et al. 2004; (9) Gonzalez et al. 2006. 
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